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ABSTRACT 

There are two main tidal effects which can act on the Local Group dwarf spheroidals 
(dSphs): tidal stripping and tidal shocking. Using N-body simulations, we show that 
, tidal stripping always leads to flat or rising projected velocity dispersions beyond a 

■ critical radius; it is ~ 5 times more likely, when averaging over all possible projection 
Q\ | angles, that the cylindrically averaged projected dispersion will rise, rather than be 
lO ■ flat. In contrast, the Local Group dSphs, as a class, show flat or falling projected 
r**"* ' velocity dispersions interior to ~ 1 kpc. This argues for tidal stripping being unim- 
portant interior to ~ 1 kpc for most of the Local Group dSphs observed so far. We 
show that tidal shocking may still be important, however, even when tidal stripping is 
not. This could explain the observed correlation for the Local Group dSphs between 
central surface brightness and distance from the nearest large galaxy. 

These results have important implications for the formation of the dSphs and for 
Q-f cosmology. As a result of the existence of cold stars at large radii in several dSphs, 

q a tidal origin for the formation of these Local Group dSphs (in which they contain 

' no dark matter) is strongly disfavoured. In the cosmological context, a naive solution 

£jj \ to the missing satellites problem is to allow only the most massive substructure dark 

& ■ matter halos around the Milky Way to form stars. It is possible for dSphs to reside 

within these halos (~ 1O 1O M0) and have their velocity dispersions lowered through 
the action of tidal shocks, but only if they have a central density core in their dark 

■ matter, rather than a cusp. A central density cusp persists even after unrealistically 
5h ' extreme tidal shocking and leads to central velocity dispersions which are too high to 

be consistent with data from the Local Group dSphs. dSphs can reside within cuspy 
dark matter halos if their halos are less massive (~ 10 9 Mq) and therefore have smaller 
central velocity dispersions initially. 
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1 INTRODUCTION 

The Local Group of galaxies provides a unique test-bed 
for galaxy formation models and cosmology. The abun- 
dance, spatial distribution and internal mass distribution, 
of satellites within the Local Group ca n provide sensitive 
tests of cosmological predictions (see e.g . |Mc^re_et_aL 1999, 
iKlvpin et ai]ll999t iKravtsov et al1l2004L iMaver et alJl2001al 
and lMaver et alJl2001blT However, central to such studies is 
an understanding of what the Local Group satellite galax- 
ies are. Observationally, they are usually split into three 
types: dSph galaxies, which typically have old stellar popu- 
lations, are spheroidal in morphology, lie close to their host 
galaxy x and are devoid of HI gas; dlrr galaxies, which have 
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younger stellar populations, irregular morphology, lie fur- 
ther away from their host galaxy and contain significant HI 
gas; and the transiti on galaxies, w hich are in between the 
dSph and dlrr types <Matedll99gh . 

Since the dSph galaxies lie, in general, closer to 
their host galaxy (there are notable exceptions - see e.g. 
iGrebel et al] l2003) it has often been argued that the tidal 
field of the host galaxy has an important ro le to play in 
the for mati on and evolut i on of t hese galaxies. IMaver et alJ 
fcOOTal) and lMaver et alJ fcOOlbl) have suggested that all of 
the Local Group satellite galaxies started out looking more 



We use the terminology 'host galaxy' throughout this paper to 
refer to either the Milky Way or M31 depending on which of these 
is closer to the satellite being discussed. 
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like the dlr rs, the dSphs then for ming through tidal trans- 
formations. iKuhn fc Milierl ]l989l) and IKroupal il997h have 
argued that the dSphs might have formed in the tidal tail of 
a previous major merger; in their mod el, dSphs con t ain no 
dark matter and are nearly unbound. IStoehr et al] J2002t) 
have argued that the dSphs inhabit the most massive sub- 
structure dark matter halos, of mass ~ 10 10 Mp), predicted 
by co smo logical N-body sim ulations (see e.g. iKlypin et al] 
1999 and Moo re et al ] |l99Sl) : in such a scenario, tides are 
unlikely to have significantly altered the visible distribution 
of stars in these galaxies. 

In this paper we compare a suite of N-body simulations 
of dSph galaxies orbiting in a tidal field, with the latest 
data from the Local Group dSphs (see Figure . In this 
way, we can constrain the importance of tidal effects for 
these galaxies and break the degeneracies between the above 
formation scenarios. 

The effects of tides on the Local Group dSphs has been 
investigated exten s ively i n the literature, b oth dynamically 
jAguilar fc White! Il986t lOh et al] ll99Sl iPiatek fc Prvorl 
[ 19951. lHavashi et al] l200l iKazantzidis et al] 120041 and 
lMashch grikci_e_LaDl20Q5|), and within a cosmological frame- 
work jKravtsovetaLl2004) . Thi s present work c o mplem ents 
that of previous authors. Like lHavashi et al.l (I2003T) and 
IKazantzidis et all (|2C)04) , we perform simulations which are 
not fully cosmologically consistent, but which have the ad- 
vantage that the simulated satellites have very high reso- 
lution. Unlike these previous authors, however, our dSph 
galaxies contain both stars and dark matter. This allows us 
to compare our results more easily with the observational 
data from the Local Group dSphs. 

We consider three possible scenarios for the tidal evo- 
lution of the Local Group dSphs. In the first (model A), 
we revisit the hypothesis that the Local Group dSphs con- 
tain no dark matter and that their high velocity dispersions 
arise instead from the action of tides 2 . In the second sce- 
nario (models B and C), we investigate the hypothesis that 
the dSphs are dark matter dominated, but that their outer 
regions have been shaped by tidal effects. In this case we ask 
whether tides could directly cause the drop in projected ve- 
locity disper sion in the outer region s of the dSphs, recently 
observed bv lWilkinson et al] i2004l) ; and see also Figure^ 
In the third scenario (model D), we investigate a scenario 
where dSphs were much more massive in the past. In this 
case, tidal shocks cause their central density and velocity 
dispersions to fall significantly over a Hubble time, becom- 
ing consistent with current data at the present epoch. It 
is important to state explicitly here that we do not wish 
to model any one specific Local Group dSph, which would 



2 A dark-matter free model for the formation of the Local Group 
dSphs is often referred to in t he litera t ure as the 'tidal model'. 
Previous tidal models, such as IKroupal ll997t) . tried to use tides 
to explain even the central velocity dispersion of the Local Group 
dSphs. As a result, in their models, after a Hubble time, the 
dSphs became fully unbound. In the models we present here, we 
use unrealistically large central mass to light ratios in order to fit 
the high central velocity dispersions of the Local Group dSphs. 
This means that tides need not work quite so hard to reproduce 
the large outer velocity dispersions. If our tidal model is r uled out, 
more extreme tidal models such as that investigated by IKroupal 
119971) are also ruled out. 



require an extensive search through parameter space and 
probably not be relevant when taken out of a cosmological 
context. Rather, we wish to investigate the generic effects of 
tides and how they reshape a dSph. 

This paper is organised as follows: in section [5] we 
briefly outline some key analytic results for tidal stripping 
and tidal shocking - the two dominant tidal effects which 
act on a satellite orbiting in a host galaxy potential. The 
analytic results and definitions therein will be referred to 
throughout this paper and are useful in both understanding 
and testing the simulation results. In section |31 we describe 
our numerical method for setting up the initial conditions 
and integrating the orbit of the satellite in the fixed potential 
of the Milky Way. In section 21 we present the results from 
a set of representative N-body simulations simulations. We 
demonstrate that tidal stripping will always lead to a flat or 
rising projected velocity dispersion. In section we discuss 
our results in the context of cosmology; we place mass limits 
on the Local Group dSphs. Finally, in section |S| we present 
our conclusions. 



2 THEORETICAL BACKGROUND 

A satellite in orbit around a host galaxy will experience 
two main tidal effects which will reshape the stellar and 
dark matter distributions: tidal stripping and tidal shock- 
ing. Since both of the mechanisms will be referred to many 
times throughout this paper, it is instructive to briefly out- 
line the salient properties of each in this section. 



2.1 Tidal stripping 

We discuss the theory o f tidal stripping in much more detail 
in a companion paper. iRead et al] i2005t) . Here we briefly 
summarise the main results of that work. Tidal stripping 
refers to the capture of stars by a host galaxy from a satel- 
lite galaxy beyond a critical radius: the tidal radius, r t . In 
general, the tidal radius depends upon four factors: the po- 
tential of the host galaxy, the potential of the satellite, the 
orbit of t he satellite and, wh ich is new to the calculation pre- 
sented in |Rea^^t^L|jj2005|K the orbit of the star within the 
satellite. In lRead et al] (|2005), we demonstrate that this last 
point is critical and suggest using three tidal radii to cover 
the range of orbits of stars within the satellite. In this way 
we show explicitly that prograde star orbits will be more eas- 
ily stripped than radial orbits; while radial orbits are more 
easily stripped than retrograde ones. 

In general, these t hree tidal radii mus t be calculated nu- 
merically as outlined in lRead et al] i2005t) . We will use these 
more accurate numerical calculations in the presentation of 
tidal radii later on in this paper. However, for the special 
case of point mass potentials, rt reduces to the following 
simple form: 
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where x v and e are the pericentre and eccentricity of the 
satellite's orbit, M s is the total mass of the satellite, M g is 
the mass of the host galaxy and a — —1, 0, 1 parameterises 
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Figure 1. Surface brightness profil es (left) and projected velocity dispersions ( right) for three Local Group d Sphs observed with goo d 
stellar kine matics. Data taken from [ W ilkinson et alJ j2004f >. fKlevna et alJ feOO^ . Iirwin fc Hatzidimitrioul <199fll) and lPalma et ah! 12003); 
but see also Martmcz-Delftado et al. (2001) for deeper surface brightness data for UMi. Notice the similarities between each of the dSphs: 
the surface brightness profiles are roughly exponential, while the projected velocity dispersions are approximately flat with a sharp fall-off 
at large radii (~ lkpc). 



the orbit of the star within the satellite. From equation 
we can see that stars on prograde orbits (a = 1) are more 
easily stripped than those on radial orbits (a — 0), which are 
more easily stripped than those on retrograde orbits (a = 
— 1). Notice that a = recovers the standard iKind dl962T) 
tidal radius. Other analytic solutions also exist for power- 
law and a restricted class of split power-law density profiles 

Over long times, orbital transformations cause a con- 
vergence of these three tidal radii on the prograde stripping 
radius and leads to the onset of tangential anisotropy be- 
yond this point. 

A final key result from iRead et al] (120051) is that the 



standard intuition that the tidal radius depends only (within 
a factor) on the densities of the host galaxy and satellite 
galaxy is only valid for stars within the satellite which are on 
radial orbits. For general star orbits, their stripping radius 
will depend upon the mass distribution within the satellite 
and not just the enclosed mass. 



2.2 Tidal shocking 

A second important effect of tides is tidal shocking. Tidal 
shocks can be caused either as a sate llite plunges through 
the disc of the host galaxy (see e.g. lOstriker et aljll972l) 
or moves on a highly eccentric orbit through the galactic 



4 Read et. al. 



Model 


p* 


PDM 


iV* 


£* (kpc) 


N DM 


£,dm (kpc) 


Orbit 


Time ( Gyrs ) 


7 


A 


P: 5, 0.23 


None 


10 5 


0.01 






9kpc,84kpc,-34.7° 


10.43 


0.0011 


B 


P: 0.072, 0.23 


P: 14, 0.5 


10 s 


0.01 


10 6 


0.01 


23kpc,85 kpc,-34.7° 


4.5 


0.00014 


C 


SP: 0.072, 0.1, 


SP: 90, 1.956, 1 


10 5 


0.01 


10 6 


0.01 


23kpc,85 kpc,-34.7° 


4.5 


0.15 


D 


SP: 0.072, 0.1, 


SP: 10 3 , 1.956, 1 


10 5 


0.01 


10 6 


0.01 


6.5kpc,80kpc,7.25° 


8.85 


0.83 



Table 1. Simulation initial conditions for models A-D. The columns from left to right show the model (labelled by a letter in order of 
discussion), the stellar density profile, p, , the dark matter density profile, pdm, the number of stars, JV*, the force softening for the 
stars, £», the number of dark matter particles, iVrjMi the force softening for the dark matter, £dmj the dwarf galaxy orbit, the output 
time in Gyrs at which we show the results from the model and the strength of tidal shocks, 7 (see section 12.21 . The density profiles 
are either (P)lummer with parameters: mass (10 7 Mq), scale length (kpc) or (S)plit (P)ower law with parameters: mass (1O 7 M ), scale 
length (kpc) and central log slope, a (see equations II ll and !13l for more details). The orbit for the dwarf galaxy is given by its pericentre, 
apocentre and inclination angle to the Milky Way disc. 
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Figure 2. Equilibrium tests for model A: The upper panels are for the stars, while the lower panels are for the dark matter. The left 
panels show the cumulative mass profile, the middle panels show the velocity anisotropy and the right panels show the fractional change 
in energy as a function of time. The blue lines are for the analytic initial conditions, the black lines are for the simulated initial conditions 
with 10° dark matter and 10 5 star particles, the red lines are for the simulated in itial conditions evolved for 5 Gyrs, while the green lines 
are for a simulation set up using the Maxwellian approximation <Hernauistfll993ft and evolved for 5 Gyrs. 



centre (see e.g. ISpitzerlll987l>. Both s c enarios are discussed 
in d etail bvlGnedin fc Ostrikerl dl997i) . lGnedhi et al.l (Il999l) 
and lGnedin et al.1 <ll99sl) ~ 

In the impulsive limit, the mean energy injected into 
the satellite at the r.m .s. radius, r, for disc shocks is given 
by (iGnedin et alJll999D : 

ro^=^Sf (l + x 2 )- (2) 



where: g m is the maximal disc force along the z direction (in 
all of the calculations presented in this paper it is a assumed 
that the plane of the Milky Way disc lies perpendicular to 
the 2-axis); V z is the 2 component of the satellite velocity at 
the point it passes through the disc; x — cut (lu is the angular 
velocity of a star within the satellite at r; t ~ 2H/V Z is the 
typical shock time for a disc of scale height, H); and a is the 
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adiab atic correction exponent which is ~ 5/2 for fast shocks 
jGnedin et aUll99Sb . 

Similarly, the mean energy injected for an eccentric or- 
bit in a spherical system (here an isot hermal potential fo r 
the host galaxy is assumed) is given bv JOnedin et alll99ifl : 



AE. 



sph 



RpVr, 



(i + x 2 y 



(3) 



where vo is the asymptotic circular speed of test particles at 
large radii in the host galaxy potential, R p is the pericentric 
radius, V p is the satellite velocity at R p and x is as above, 
but the shock time is now given by, r ~ ivRp/Vp. 

We make the somewhat crude assumption that the total 
energy injected into the satellite due to a single shock at f 
is then given by the sum of the above two terms. 

Assuming that the satellite is in virial equilibrium be- 
fore the tidal shock and that it has an isotropic velocity 
distribution, the total specific energy of the satellite, E, is 
given by: 



E 



-T ■ 



1 



— v 2 (r) 

2 V ' 2p s {r) 



ps{r) ^dr 



(4) 



where T is the satellite's total specific kinetic energy, 
v 2 (r) is the velocity dispersion, p 3 (r) and M„(r) are the 
mass and density profile of the satellit e and the right hand 
side follows from the Jeans equations (Bin nev fc Tremainel 
Il987h . 

We assume that tidal shocks are important for the satel- 
lite if the energy injected is comparable to the total energy 
of the satellite initially. This gives: 



A-Esph + Ai3di S c 



v 2 (r) 



(5) 



where 7 parameterises the importance of tidal shocks: tidal 
shocks are important for 7 ~ 1 and unimportant for 7< 1. 
Notice from equations |5] 01 an d |S] that the effectiveness of 
tidal shocks depends on the mass density of the satellite, the 
potential of the host galaxy and on the orbit of the satellite. 
The strength of such shocks is very sensitive to the value 
of R p . For the Milky Way potential employed in this paper, 
both disc shocks and spherical tidal shocks are negligible for 
orbits which never come closer than ~ 20kpc. 

The value of 7 for each of the four models presented in 
this paper (labelled A-D) are given in Tabled Model D is a 
case of particular interest. For this model, 7 ~ 1 and so tidal 
shocks are very important. Yet the tidal radius (as obtained 
using equation lies at the edge of the light distribution. 
Thus it is possible that some of the Local Group dSphs could 
be in a regime where their visible light is strongly affected 
by tidal shocks but only very weakly perturbed by tidal 
stripping. 



3 THE NUMERICAL TECHNIQUE 
3.1 Setting up the initial conditions 

3.1.1 The dSph galaxy 

The initial conditions for the dSph gala x y we re set up 
in a similar fashion to IKazantzidis et all i2004h . but for 
two-component spherical galaxies comprising a dark matter 



halo and some stars. Since this procedure is non-trivial we 
briefly outline the major steps involved here: 

(i) The particle positions for each component are re- 
alised using standard accept/reject techniques with 
numerically cal culated compariso n functions and analytic 
density profiles JPress et al.lll992l) . 

(ii) We us e isotropic, spherica l ly sy mmetric, distribution 
functions feinnev fc Tremainel Il987l) . For each separate 
component, i, with density, pi, in the total relative 3 grav- 
itational potential of all components, ip, the distribution 
function is given by the Eddington formula: 



m = 



x/8tt 2 



d 2 pi dip 



dtp 2 y/e — ip 



1 f dpi 
^fe\dlp 



1/1=0 



(6) 



where e is the specific energy. Note that, for any choice of pi 
and ip which are well behaved as r — *■ 00, the second term 
in th e brackets in equation [5] vanishes feinnev fc Tremainel 
I1987D . 

For a general two component system /(e) must be cal- 
culated numerically. In lKazantzidis et al.H2004h . where only 
single component systems (dark matter only) were consid- 
ered, the term could be calculated analytically 4 . Here, 
where in general ip = tpi + ip2 + ... + ipn for an n component 
system, this must be calculated numerically. We consider 
the special case where ipi and pi are known analytically for 
each component. In this case we can write: 



d 2 pi 
dip 2 



d 2 pi 
dr 2 



dr dpi d I dr 
dip dr dr \ dip 



(7) 



d 

This means that we can calculate in general 

analytically for a given value of r. This avoids the need 
for a noisy double numerical differential. However, we do 
still need to generate a look-up table to solve the inversion 
for r(ip). Once this is done, however, we may perform 
the integral in equation [§] and c reate a log look-up tabl e 
for /(e) over a wide range of e as in lKazantzidis et al.N2004) . 

(iii) Once the distribution function, /(e), is known for 
each component, we can set up p article velocities using 
standard accept/reject techniques flPress et al.lll99^1 . We 
take advantage of the fact that the maximum value of the 
distribution function at a given point will be given by the 
maxi mum value of e; that is for v = l|Kuiiken fc Dubinskil 
Il995h . 



This may seem like an unnecessarily large amount of 
effort to go to, when simple numerical techniqu es ahead; 
exist for setting up multi-component g alaxies ( Hcrnquis 
1993). However, as pointed out by IKazantzidis et al 
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2004), using approximate methods such as the Maxwellian 
approximation will lead to a small evolution away from 
the initial conditions. This can lead to significant radial 
velocity anisotropies being introduced at large radii. In a 



3 The relati ve gravitation potent i al, il> = 
as denned in lBinnev &: Tremainel Il987l) ■ 

4 In fact, even for single component models, 
for a few special cases. 



-$ + $01 used here is 
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is only analytic 
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Figure 3. Initial conditions for models A-D. The blue data points are for the Draco dSph galaxy taken from Iwilkinson et alJ (2004). 
The black, red, green, purple lines are for models A,B,C,D respectively. Recall that in model A we have no dark matter, hence there is no 
black line in the right panel plot of the dark matter distribution. 



companion paper dRead et al.1 120051) we show that such 
anisotropies lead to significantly increased tidal stripping. 
This is undesirable when we wish to perform controlled 
numerical simulations of tidal stripping. 

In Figure H we show an equilibrium test for our two- 
component model B (the other models showed similar re- 
sults). The upper panels are for the stars, while the lower 
panels are for the dark matter. The left panels show the cu- 
mulative mass profile, the middle panels show the velocity 
anisotropy and the right panels show the fractional change 
in energy as a function of time. We use a modified anisotropy 
parameter, /3*, given by: 

* = (8) 

where v 2 and v 2 are the radial and tangential velocity dis- 
persions respectively. This definition gives /?* = — 1 for pure 
circular orbits, /3* = for isotropic orbits, and /?* = 1 for 
pure radial orbits 5 . 

The blue lines in Figure [5] are for the analytic initial 
conditions, the black lines are for the simulated initial con- 
ditions with 10 6 dark matter and 10 5 star particles, the red 
lines are for the simulated initial conditions evolved in iso- 
lation (i.e. with no external tidal field) for 5 Gyrs, while the 
green lines are for a simulation set up using the Maxwellian 
approximation jHernauisdli99il and evolved for 5 Gyrs. 

The most important point to t ake away from Figu rcElis 
that we can confirm the findings of Kaz antzidis et alJ (2004) 
for our two-component model. In the Maxwellian approxi- 
mation the dark matter evolves away rapidly and signifi- 
cantly from the initial conditions. In particular, significant 
radial anisotropy is introduced to the dark matter at large 

5 The standard anisotropy parameter, /3, is usually defined as 
P = l-vj/v 2 iBinnev fc Tremainell9'87h . This definition is useful 
for obtaining analytic solutions to the Jeans equations. However, 
it can be readily seen that for circular orbits j3 — > — oo. In this 
paper it is useful to have a definition of anisotropy which is well- 
behaved for all orbits and which is symmetrical around isotropic 
orbits. 



radii (see green line, bottom right panel). In contrast, the 
distribution function generated initial conditions are very 
stable over the whole simulation time. 

Notice that the energy both in the stars and dark matter 
is conserved to better than 2% over the whole simulation 
time. 

3.1.2 The host galaxy potential 

We used a host galaxy pote ntial chosen to provide a 
good fit to the Milky Way llLaw et al.l 120051) . with a 
Miyamoto-Nagai potential for the Milky Way disc and bulge 
jNaeai fc Mivamotclll97ri . and a logarithmic potential for 
the Mi lky Way dark matter halo. These are given by respec- 
tively feinnev fc Tremainelll98'it) : 

(R,z)= ~ GMd (9) 

^/R 2 + {a + ^/z 2 + b 2 ) 2 

where M — 5 x 1O 1O M0 is the disc mass, a = 4kpc is the 
disc scale length and b — 0.5 kpc is the disc scale height, 
and: 

'I'logM = t^o m {Rc + T ' 2 ) + constant (10) 

where R c — 4.1 kpc is the halo scale length and vo = 
220km/s is the asymptotic value of the circular speed of 
test particles at large radii in the halo. 

Note that with the choice b = 0.5 kpc, we do not con- 
sider the case of maximal disk shocking which would occur 
for b ^ (see equation 0. However, given that the amount 
of disc shocking is also degenerate with the orbit of the dSph 
which is poorly constrained, our simulations span a range 
from typical to extreme levels of shocking. In addition, none 
of our conclusions would be altered by stronger disk shocks. 

3.2 Force integration and analysis 

The initial conditions we re evolved using a v ersion of the 
GADGET N-body code JSpringel et al] 1200 ill modified to 
include a fixed potential to model the host galaxy. 
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A wide range of code tests were performed using dif- 
ferent softening criteria, force resolution and timestep cri- 
teria and the results were found to be in excellent agree- 
ment with each other. We also explicitly ch ecked the code 
against output f rom two other N-body codes : ll5ehnerj i200(J) 
and NBODY6 llAarsetblll999l) . and found excellent agree- 
ment. Our f o rce so ftening was chosen using the criteria of 
IPower et alJ (120031) and is shown in Tabled For the equi- 
librium tests, our simulations were found to conserve total 
energy to better than one part in 10 3 over the whole simu- 
lation time of 5 Gyrs. 

Finally, it is not trivial to mass and momentum centre 
the satellite when performing analysis of the numerical data. 
An incorrect mas s centre can lead to spurious density and 
velocity features jPontzen et al.ll2005li . We use the method 
of shrinking sp heres to find the m ass and momentum centre 
of the satellite JPower et alJl2003h . 



3.3 The choice of initial conditions: models A-D 

We chose four models for our initial conditions labelled A-D 
as shown in Tableland Figure |3] In all models the density 
distribution for the st ars and dark matter were given either 
by Plummer spheres feinney fc Tremainej|l98 7|) or by Split 
Powe r la w profiles (S P) (lHernauistll99dlSahall99llbehnenl 
Il993l and IZhaclll99rJl . The density-potential pairs for these 
distributions are given by respectively: 
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where M and a are the mass and scale length in both cases, 
and a is the central log-slope for the SP profile (note that 
the SP profile always goes as psp oc r -4 for r ^> a). 

As detailed in section^ each of these models was chosen 
to test a different scenario for the evolution of dSphs in the 
presence of tides. In all of the models, the dSph was placed 
on a plunging orbit which took it through the Milky Way 
disc. This ensured that both tidal stripping and shocking 
were important (see section The initial conditions, orbit 
and output time for each of the models are given in Table 
Models B and C placed the dSph on an orbit which is 
consisten t with current cons traints from the proper motion 
of Draco llKlevna et alJ200lTl . However, these proper motion 
measure ments are notoriou sly difficult and the errors large 
(see e.g. iPiatek et all 12002ft . Thus in models A and D, we 
also consider more extreme (i.e. more radial) orbits. 



4 RESULTS 

The results for models A-D are presented in Figures 2] and 
|^| In Figure^] we show the projected surface brightness (left 



panel), projected velocity (middle panel) and projected ve- 
locity dispersion 6 (right panel) for the stars. For the pro- 
jected velocity dispersion, the errors are determined from 
Poisson statistics. Models A-D are presented down the page. 
In Figure we show the projected surface brightness profile 
for the stars (left panel), the true velocity dispersion of the 
stars (middle panel) and the density profile of the dark mat- 
ter (right panel). For clarity on the plots, we compare the 
results of our simulations to data from the Draco dSph only, 
on the assumption that Draco is representative of the most 
dark matter dominated Local Group dSphs (see section 
and Figure . In Figure |H] we show the effect of different 
viewing angles in projection on the sky, using model A as 
an example. In Figure|H|we show the mass loss, due to tides, 
as a function of time for all four models A-D. 



4.1 The no-dark matter hypothesis: model A 

In model A, we test the hypothesis that we can eliminate the 
need for dark matter in dSph galaxies by instead employing 
tidal heating to raise the velocity dispersions. One of the 
appealing aspects of such a model is that it produces tidally 
distorted isophotes, such as those which may have been seen 
in some dSphs (s ee e.g. iMartmez-Delgado et al.ll200ll and 
iPalma et aT]l2003l) . 

In Figure 2] top left panel, we can see that distortions 
occur naturally in such a tidal model and the S-shaped 
isophotes caused by tides can be clearly seen in this pro- 
jection. However, once the kinematic data are taken into 
account, it is clear that such a model runs into serious 
difficulties. The projected velocity shows a strong gradient 
(10 km/s/kpc) due to the near-dominant tidal arms, whereas 
no such velocity grad ients have been observ ed in the Local 
Group dSphs so far JWilkinson et al.ll2004l) . Furthermore, 
the projected velocity dispersion rises out towards the edge 
of the light, whereas the Local Group dSphs like Draco 
have flat or falling projected velocity dispersions (see e.g. 
IWilkinson et al.l2004l.lMunoz et all2005llKlevna et al.l2004l 
and Figure 0. It is possible, however, that such tidal fea- 
tures could be hidden in projection. Figure|S]shows the effect 
of different viewing projections on the sky for model A. The 
top and bottom panels show two different projections: the 
top is a more typical projection, while the bottom, which 
kinematically hides one of the tidal tails, is more rare. By 
rare, we mean that such obscurations only occur for ~ 1/5 
of all projection angles. While the orbit of any given dSph 
would rule out some projection angles, when studying the 
generic effect of tides on all of the Local Group dSphs, it 
seems reasonable to consider such random projections. 

Notice from the middle left panels, that the primary 
cause of the rising velocity dispersion in projection is the 
cylindrical averaging which picks up both tidal tails. This 
is why, when a suitable projection is found which hides one 



6 Note that throughout this paper, all velocity dispersions are 
calculated using the local mean velocity; this is in contrast to ob- 
servational techniques which often use the global mean for all the 
stars. The use of the global mean is only justified where there are 
no velocity gradients measured; this is the case for the observa- 
tions. It is not the case, however, for some of the models presented 
here where velocity gradients are induced through tides. 
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Figure 4. Projected surface density (left), velocity (middle) and velocity dispersion (right) for models A (top) through D (bottom). The 
grey contours show particle numbers per unit area (left) and velo cities in km/s (middle). Over-plotted on the projected velocity dispersion 
plot (right) are the data for the Draco dSph galaxy taken from Iwilkinson et al ] l2004h . The blue boxes overlaid on the contours mark 
approximate current observational limits. The projections shown are typical; different projections for model A are shown in Figure Itjl 



of the tails, the velocity dispersion becomes fiat rather than 
rising (bottom panels). In this case the velocity gradients are 
visible (middle right panel) only at large radii. Interior to 
~ 1 kpc, the velocity gradient is less than ~ 5 km/s. Such a 
small vel ocity gradient would b e difficult to detect observa- 
tionally iWilkinson et al.ll2004l) . However, even this special 
projection is inconsistent with the data from Draco, UMi 



and Sextans which all show falling velocity dispersions in 
projection (see Figure 0. 

It is worth recalling at this point that we are not try- 
ing to explicitly model the Draco, UMi or Sextans dSphs. 
Instead we show that the generic effect of tides is to pro- 
duce velocity gradients and flat or rising projected velocity 
dispersions. That this is not seen in Draco, Sextans or UMi 
suggests that while tidal stripping may be important be- 
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Figure 5. Projected surface brightness profile of the stars (left panel), velocity dispersion of the stars (middle panel) and the de nsity 
profile of the dark matter (right panel). Over-plotted are three theoretical tidal radii calculated as described in lRead et alj 120051) : the 
blue/green/red vertical lines are for the prograde/radial/retrograde tidal radii respectively. The black lines show the initial conditions, 
while the red lines show the evolved profiles at the output times given in Tabic In the velocity dispersion plot, the solid, dotted and 
dashed lines show the r, 9 and <j> components of the velocity dispersion respectively. Models A-D are shown from top to bottom. Model 
A did not contain any dark matter which is why no plot is shown for the dark matter in this case. The tidal radius for prograde star 
orbits (vertical blue line) cannot be seen for model A since it lies at very small radii (9pc). 



10 Read et. al. 








O.S 1.0 l.S 2.0 J,S iJB 

ft {"pO 



15 

1 



c 10 

'a 




6.S 1.0 l.S 2.0 J,S -*jO 

ft {"pO 



Figure 6. The effect of viewing projection for model A. Top panels show a typical projection; bottom panels show a rare projection 
which kinematically hides one of the tidal tails. From left to right panels show: projected surface density, individual projected star 
velocities, mean projected velocity and projected velocity dispersion. The grey contours show particle numbers per unit area as in Figure 
PI Over-plotted on the projected velocity dispersion plot (right) are the data for the Draco dSph galaxy taken from Wilk inson et al] 
12004). The blue boxes overlaid on the contours mark approximate current observational limits. Note that the velocity gradients shown 
here are averaged over the y and z coordinates and are necessarily different to the two dimensional gradients shown in Figure HI 



yond the edge of the currently observed light (~1 kpc), it is 
unlikely to have been important within this radius. 

An interesting feature of model A is that the small 
pericentre of the satellite (9 kpc) causes strong tidal shocks 
which over 10 Gyr lower the central density of the stars by 
a factor of ~ 2. This can be seen in Figure El top p anels. 
These tidal shocks, as discussed in lRead et alJ (120051) wash 
out the tidal features which would otherwise be expected at 
the analytic tidal radii (blue, red and green vertical lines in 
Figure UJ. However, the effects of tides can still be seen in 
the velocity dispersions. Notice that tangential anisotropy, 
which should be present beyond the prograde stripping ra- 
dius is present at all radii over the satellite (see section l2~T]l . 
For model A, the prograde stripping radius lies at just 9pc 
and hence does not show up in Figure Q3 



4.2 Weak tides: models B and C 

In models B and C we consider the effect of weak tides on 
the Local Group dSph galaxies. In this case we imagine that 
tides gently shape the outer regions of dSphs over a Hubble 
time. 

In Figure [I] second and third panels from top, we can 
see that tides have little affected models B and C. The sur- 
face brightness distributions are nearly spherical, as in the 
initial conditions, and there are no visible velocity gradients 
across the galaxy. Tidal heating of the outermost stars is 
present in both models from about ~ 2 kpc outwards. No- 
tice that in neither model is there a drop in the projected 



velocity dispersion. The decline in model A was present in 
the initial conditions and is not as sharp as that seen in the 
data from Draco. 

This can also be seen in Figure^] second and third pan- 
els from top. Notice that, as in model A, tangential velocity 
anisotropy appears at the prograde stripping radius (vertical 
blue line) in both models B and C. Notice also the action of 
tidal shocks, particularly on the dark matter in model C. As 
the satellite is puffed up by shocks, its velocity dispersion 
at all radii is lowered. This may seem counter-intuitive at 
first: the satellite is heated by tidal shocks and yet its veloc- 
ity dispersion falls. However, it is the subsequent expansion 
of the satellite after the injection of energy from the tidal 
shock which causes a drop in its velocity dispersion at a 
given radius . 

Notice further that it is only in projection that the ve- 
locity dispersions appear hot beyond the tidal radius. The 
true velocity dispersions are still gently falling, as in the ini- 
tial conditions, at this point. This occurs for two reasons. 
First, and this is usually the dominant effect for a given 
projection, the cylindrically averaged tidal tails cause the 
projected velocity dispersion to rise (see Figure ^J. Second, 
the onset of tangential anisotropy makes the outermost stars 



7 In projection this becomes a little more complicated since hot 
stars which have moved out to large radii can inflate the projected 
velocity dispersion. This docs not happen in practice, however, be- 
cause these stars move beyond the tidal stripping radius, become 
unbound and move away from the satellite. 
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Figure 7. Histogram of the stellar velocities in a 0.2 kpc bin at lkpc for model B (left) and model C (right). Over-plotted are Gaussian 
fits in both cases (red lines). 



appear hot in projection. This second effect is only impor- 
tant when the tidal tails are viewed along the line of sight, 
thereby minimising the extent to which they can be observed 
in projection. 

This is interesting: in all models tidal stripping leads 
to a rising or flat projected velocity dispersion. There is 
never a sharp decline such as that observed in the Draco, 
UMi and Sextans dSphs (see section . We now discuss 
why tidal stripping always leads to hot stars in projection 
by considering three mechanisms by which tidal stripping 
might have led instead to the opposite - namely a cold 
outer point in the projected velocity dispersion. 

(i) The outer-most point is the result of a chance pro- 
jection effect between the galaxy and its tidal tail along the 
line of sight. 

We explicitly checked for this by searching over 5 de- 
gree projections every half orbital period for both models B 
and C and found no sharp features at any output time. If a 
projection effect could cause such a sharp drop it would be 
very rare and, while one could envisage this for one Local 
Group dSph galaxy, the fact that such a sharp feature 
has been observed in at least three renders this possibility 
unlikely. 

(ii) The dSph galaxy is on a significantly plunging or- 
bit. The outermost stars are heated up by a disc shock as 
the dSph plunges through the Milky Way disc. The dSph 
galaxy then moves out towards apocentre and the hot outer 
stars escape leaving a cold population behind. 



This hypothesis was presented in I Wilkinson et alJ (12004 
and formed one of the motivations for this current work. 
The output time chosen for models B and C has been 



specifically selected to test this hypothesis. We choose a 
time when the dwarf galaxy has spent a large amount of 
time away from the Milky Way disc and the outermost stars 
heated by the last passage have had time to escape. The 
hypothesis is falsified for three reasons. First, for the case 
of models B and C, tidal shocking replenishes stars beyond 
the tidal radius and so there is always a source of hot 
tidal stars which mask any cold populations. Secondly, the 
image of a sharp edge to the satellite beyond which a cold 
population of star s might exist is no t correct. As shown in a 
companion paper dRead et alJl2005t and see also section 
the tidal radius of the satellite depends on the orbit of the 
stars within the satellite. This leads to a continuum of tidal 
radii and, therefore, no sharp edge to the satellite. Finally, 
even if there were such a sharp edge, and tidal shocking did 
not refill stars beyond the tidal radius, there would be on 
average more stars on circular orbits than radial orbits at 
the tidal radius and the distribution would still appear hot 
in projection. We explicitly checked that this is indeed the 
case by considering some more extreme orbits with very 
large apocentres (260 kpc), and pericentres large enough to 
make the effect of tidal shocking negligible (35 kpc). Even 
in these cases the projected velocity dispersions are flat or 
rising, for the reasons presented above. As a final check, we 
considered also only the bound stars. Plotting only these 
removes the hot tidal tail, but still does not lead to cold 
outer stars in projection. 

We conclude that if such cold outer populations in 
dSph galaxies are real then they cannot have formed as 
a result of tidal effects which always produce stars which 
appear hot in projection. 

(Hi) The velocity distribution of the outer-most stars 
is highly non-Gaussian, with a significant tail of high 
velocity stars. This could lead to a large fraction of stars 
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being deemed non-members of the dSph galaxy and excluded 
from the data analysis, leading to a spurious detection of a 
cold outer point. 

To test this hypothesis, we plot the distribution of 
line of sight velocities in a 0.2kpc-wide bin at a projected 
radius of lkpc for models B and C (Figure \Q. Over-plotted 
are Gaussian fits to the binned data in both cases (red 
lines). The distributions are close to Gaussian: neither 
shows a strong asymmetric tail. Similar results were found 
at other radii. We conclude that it is unlikely that a 
significant fraction of high velocity stars would be excluded 
as members of Draco from an observational data analysis. 

4.3 Strong tidal shocking: model D 

In model D, we consider the hypothesis that the progenitors 
of dSph galaxies were much more massive than they appear 
at the present time. As can be seen in models A, B and C, 
the action of tidal shocks puffs up the satellite and lowers 
the velocity dispersion at all radii. In model D, we start with 
a satellite which has a very large central velocity dispersion 
which is not consistent with current data from the Local 
Group dSphs. We consider the case where this velocity dis- 
persion is lowered over time by the action of tidal shocks 
such that it becomes consistent with current data now. 

As we showed in section 12.21 it is possible to find a 
regime where tidal shocks dominate over stripping interior 
to the stars. The existence of such a regime is demonstrated 
numerically by model D. 

In Figure^] bottom panels, we can see that the satellite 
has been little affected by tidal stripping inside ~ I kpc. It 
retains smooth spherical contours and shows no strong ve- 
locity gradients. Beyond ~ 1 kpc, tidal stripping becomes 
significant and the outer stars appear hot in projection. 
However, in Figure |S] we can see that tidal shocking has 
been much more important than the stripping over a Hub- 
ble time. The central density of the dark matter has lowered 
by a factor of ~ 10, and similarly for the stars. The velocity 
dispersion of the stars is lower at all radii, but notice that 
the central value is still high - similar to the initial condi- 
tions. This reflects the presence initially of a central dark 
matter cusp. This allows the density at the centre of the 
satellite to become high enough that tidal shocking is less 
effective, even over very long timescales and for this extreme 
orbit (recall that the satellite pericentre for model D is just 
6 kpc). Finally, notice that despite such strong tidal shocks, 
the central dark matter cusp persists, albeit at a lower nor- 
malisation. This agrees w ith findings from previous authors 
( Kaz antzidis et al.ll20 04). It is this central dark matter cusp 
which leaves the satellite with a central velocity dispersion 
which is too large to be consistent with data from any of 
the Local Group dSphs observed so far. In a similar model 
which we ran with a central dark matter core instead of a 
cusp, much better agreement was obtained. We will return 
to this issue in section [K] 

In model D, as in models A and C, the satellite loses sig- 
nificant mass over a Hubble time as a result of tidal stripping 
and shocking. This is shown in Figure [S] where we compare 
the bound mass as a function of time for each of the mod- 
els A-D. At the output time shown, the bound mass of the 
satellite in model D has been reduced by a factor of ~ 10. 
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Figure 8. Total bound mass (stars and dark matter) as a function 
of time for models A-D. 



This explains why the final velocity dispersion matches the 
data from Draco, while the initial dispersions were much 
too high. Such mass loss continues steadily over the whole 
simulation time due to the action of shocks. 



5 DISCUSSION 

5.1 Sagittarius: a tidally stripped dSph galaxy 

We know of one Local Group dSph which is definitely u nder- 
goin g tidal stripping: the Sagittarius dwarf lllbata et al.l200ll 
and lLaw et al]l2005tl . Its tidal streams ar e quite cold, giv - 
ing a flat velocity dispersion of ~ lOkm /s tLaw_ et al]l2005l) , 
consistent with its central value lllbata et~aL^j97l) . Further- 
more the velocity gradient s across its mino r axis (±lkpc) 



see m to be onl y ~ 3k m/s ll^^^^^al 



lOlszewskil <1998ft and 



Ibata et al 



119971)* 
' lll997ft 



have argued 

that these data for the Sagittarius dwarf must imply large 
amounts of dark matter, similar to the other dSphs, such 
that the effects of tides become hard to detect within the 
body of the Sagittarius dSph. An alternative view is that 
even strong tides do not produce measurable effects in the 
kinematics, even out to radii which begin to sample the tidal 
tails. As we have shown in this paper, this second view is pos- 
sible if Sagittarius is favourably aligned (see Figure [SJ. One 
might speculate that the other dSphs could be undergoing 
tidal disruption but the signs of it are masked in projection. 
Such favourable alignments seem unlikely to be the case for 
all of the Local Group dSphs. Furthermore, in at least three 
dSphs (UMi, Draco and Sextans), the velocity dispersion is 



8 The major axis results are not available in the published litera- 
ture. It is possible that large velocity gradients may be discovered 
in future along the major axis, consistent with strong tides. 
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actually falling. This is inconsistent with strong tidal effects 
interior to ~ lkpc. 

5.2 The tidal model 



lKroupa| lll997h found that tidal models for the formation of 
dSphs without dark matter can reproduce the spatial and 
kinematic observations. Central to such models is the idea 
that the dSph galaxy starts out as a high density star clus- 
ter. Through tidal shocking and stripping the star cluster is 
heated up and becomes almost completely unbound. In such 
a scenario, high projected velocity dispersions are to be ex- 
pected (c.f. model A). However, these do not correspond to a 
large mass for the satellite since the tidally heated stars are 
not bound and the Jeans equations or distribution function 
modelling do not apply. 

In ou r simulat i ons, w e find sim i lar re sults to those pre- 
sented in iKroupal l)l997l) . iKroupal 1)19971) also finds rising 
(cylindrically averaged) projected velocity dispersions for 
tidally stripped and shocked dSphs (see their Figure 10). It 
is the new d ata from the Lo cal Group dSphs, which were not 
available to IKroupal (|l997h , that rule out these dark matter 
free tidal models, not an improvement in the simulations. 
Even using favourable projections (as in Figure it is not 
possible to have a, falling projected velocity dispersion in the 
tidal model - a result which concords with other studies of 
the tidal model too ( see e. g. iFleck fc Kuhnll2003fi . 

iPiatek fc Prvorl lll995t) also argued against a dark mat- 
ter free tidal origin for the dSphs. They found significant 
velocity gradients in their sim ulations which were inconsis- 
tent with observations. As in IKroupal l)l997h . we find that 
it is possible to hide such velocity gradients using suitable 
projections. Thus it is the existence of cold outer popula- 
tions in projection which provide the hardest challenge for 
such tidal models, rathe r than the velocity gr adients alone. 

Finally, for Draco, iKlessen et alJ l)2003|) showed that 
the tidal model will not work since Draco cannot be very 
extended along the line of sight (something which is re- 
quired in any projection which disguises the velocity gradi- 
ents and significant tidal tails, but in which the velocity dis- 
persion is inflated due the inclusion (in projection) of stars 
from the tidal tail s ). Thr ough detailed modelling of Draco, 
iMashchenko et all (120051) also find that the tidal model is 
not favourable. Here, we confirm these results, and extend 
the list of tidal stripping-free dwarfs to include UMi and 
Sextans. 

5.3 Mass bounds from tides 

It is interesting to now turn the problem of the cold outer 
stars in Draco, Sextans and UMi on its head. The presence 
of such stars tells us that either the tidal radius must lie be- 
yond this point, or that the cold point was formed recently 
and has not had time to be disrupted by tides. The latter 
seems unlikely: Draco and UMi have very old stellar popula- 
tions ( £ 10 Gyrs), which suggests that any baryonic process 
which might cause such a cold point must have occurred long 
ago. As such, in this section, we discuss the implications of 
these Local Group dSphs having tidal radii ^ 0.8 kpc, which 
is the current observational edge for kinematic observations 
in these galaxies. For brevity from here on, we use the term 
'Local Group dSphs' to mean Draco, UMi and Sextans. 
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Figure 9. A lower mass bound on the Local Group dSphs as 
a function of their orbital pericentres. The red line is for highly 
eccentric dSph orbits (e = 1) while the black is for circular orbits 
(e = 0). The mass limits were obtained from the prograde tidal 
radius (see equation with a = 1). Also marked is the approxi- 
mate lower bound on the mass obtained from stellar kinematics 
and dis tribution function models for Draco (black dashed line; 
and see iKlevna et all200ll) . UMi and Sextans have similar lower 
bounds. 



From Figures [I] and [3] we can see that the tidal heating 
of stars starts at the prograde stripping radius (this is where 
the tangential velocit y anisotropy begins - see section 12.11 
and iRead et alJboOEft . Thus from here on, by 'tidal radius' 
we refer to the prograde stripping radius. 

The tidal radius is a function of the satellite and Milky 
Way potentials, and the orbit of the satellite. However, as- 
suming point masses leads to, a t most, errors of t he order of 
2-3 in the prograde tidal radius jRead et alJl2005h . Thus, re- 
arranging equationQshould give us a reasonable lower bound 
on the mass of Local Group dSphs within the tidal radius, 
as a function of their orbital pericentre. This is plotted in 
Figure |5] It is a lower bound because the tidal radius must 
be greater than the outermost measured kinematic point 
(in order to preserve a cold outer point) in the Local Group 
dSphs, which we take to be ~ 0.8 kpc. 

Since we have shown that tidal stripping is not likely to 
act interior to ~ 0.8 kpc in at least Draco, UMi and Sextans, 
this suggests that these galaxies are close to equilibrium and 
a mass bound from the Jeans equations or distribution func- 
tion modelling is sound. What, then, can we learn from a 
mass bound from tides which is degenerate with the orbit 
of the dSph and the potential of the Milky Way, both of 
which are poorly constrained? From Figure [5] we can see 
that, if the orbital pericentre of the Local Group dSphs is 
^ 30 kpc, then we learn nothing new from the tidal mass 
bound: tidal stripping and shocking are both negligible for 
all satellite masses Si 10 8 Mq - the mass bound obtained 
from distribution function modelling. However, if the peri- 
centre of the orbit is 30kpc, then the satellite must have 
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been more massive in the past than it is at present in or- 
der to have survived intact on its current orbit. Thus for 
small pericentre orbits, the tidal mass bound is a bound on 
the dSph's initial mass, before tidal stripping and shocking. 
That it is now less massive than this then means that tidal 
shocking must have acted to lower the central mass and ve- 
locity dispersion over a Hubble time (as in models C and 
D). 

5.4 Can dSphs reside in the most massive 
substructure halos? 

The dark matter halo in model D had a total mass of 
1O 1O M0, consistent with the most massive dark matter halos 
predi cted by cosmologica l simulations to surround the Milky 
Way jKlvpin et al The initial mass within 0.8 kpc for 

model D was 1O 9 M0, while after significant tidal shocking on 
its extreme orbit (with pericentre of 6.5 kpc), the final mass 
within 0.8 kpc was ~ 1O 8 M0. This final mass is consistent 
with current mass estimates for Draco, UMi and Sextans, 
but its concentration is not. The survival of the central dark 
matter density cusp in model D ensured that the central ve- 
locity dispersion of the surviving dSph galaxy was too large, 
even after extreme tidal shocking. 

Interestingly, there may be tentative evidence that the 
UMi dSph has a cen trally cored dark ma tter density profile, 
rather than a cusp jKlevna et alj feoOS). If this is the case 
then it is possible that UMi started out with total mass 
1O 1O M0 and had its velocity dispersion lowered through the 
action of tidal shocks. With a central dark matter density 
core, rather than a cusp, such shocks can lower even the cen- 
tral velocity dispersion to lOkm/s in accord with the mea- 
sured values for the Local Group dSphs. Thus it is only pos- 
sible for the dSphs to inhabit the most massive substructure 
halos if they have density profiles which are not consistent 
with simple cosmological predictions. 

A natural consequence of such strong tidal shocking is 
that we can expect a correlation between the central sur- 
face brightness of the Local Group dSphs and their peri- 
centre: dSphs which pass closer to their host galaxies will 
show lower central surface brightnesses. While the orbits 
of the Local Group dSphs remain poorly constrained and 
so their pericentres are not well known, a correlation be- 
tween central surface brightness and distance to the nea rest 
host galaxy has been obse rved jBellazzini et aljll99rJ and 
IMcconnachie fc Irwinll2C)05l) . While there are many possible 
explanations for such a correlation (for example selection 
effects), a plausible explanation would be tidal shocking. 

Finally, tidal shocks could provide the necessary physics 
to remove angular momentum from the dSphs and transform 
them from an initi ally disc-like populat ion to a spheroidal 
populat ion as in IMaver et al] ]2001all and IMaver et al] 
(2001b). We argue in this paper that the Local Group dSphs 
have not been significantly stripped interior to 1 kpc, but this 
ne ed not invalidate the tidal transformation model proposed 
bv lMaver et"aH feOOlal) . 



6 CONCLUSIONS 

We have compared the results from a suite of N-body simula- 
tions of the tidal stripping and shocking of two-component, 



spherical, dwarf galaxies in orbit around a Milky Way like 
galaxy, with recently obtained kinematic data for the Local 
Group dwarf spheroidals (dSphs). 

Our main findings are summarised below: 

(i) Tidal stripping always leads to flat or rising pro- 
jected velocity dispersions beyond a critical radius; it is 
~ 5 times more likely, over random projections, that the 
cylindrically averaged projected dispersion will rise than be 
flat. 

(ii) In several of the dSphs observed so far, there ap- 
pears to be a sharp fall-off in the projected velocity 
dispersion at large radii. If such a feature is real and not a 
statistical artifact, it would be rapidly destroyed by tidal 
stripping - tidal stripping cannot have acted interior to the 
currently visible light in these galaxies (interior to ~ 1 kpc) . 

(iii) By contrast, there exists a regime in which tidal 
shocking can be important interior to the visible light, even 
when tidal stripping is not. This could explain the observed 
correlation for the Local Group dSphs between central 
surface bright ness and distan ce fr om the nearest large 
galax y dBellazzini et al] Il996l and IMcconnachie fc Irwinl 

(iv) It is possible for dSphs to reside within the most 
massive substructure dark matter halos (~ 1O 1O M0) and 
have their velocity dispersions lowered through the action 
of tidal shocks, but only if they have a central density 
core in their dark matter, rather than a cusp. A central 
density cusp persists even after unrealistically extreme tidal 
shocking and leads to central velocity dispersions which are 
too high to be consistent with data from the Local Group 
dSphs. dSphs can reside within cuspy dark matter halos if 
their halos are less massive (~ 1O 9 M0) and therefore have 
smaller central velocity dispersions initially. 

(v) A tidal origin for the formation of the Local Group 
dSphs in which they contain no dark matter is strongly 
disfavoured. This is because such a formation scenario 
requires very strong tidal stripping, the effects of which are 
at odds with the latest data from the Local Group dSphs. 
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